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Abstract 

Bulk viscosity suppresses elliptic flow Vz, as does shear viscosity. It can thus not be neglected 
when extracting the shear viscosity from elliptic flow data. We here explore uncertainties in 
the bulk viscous contribution to viscous V2 suppression that arise from presently uncontrolled 
uncertainties in the initial value of the bulk viscous pressure and its microscopic relaxation time. 



1. Introduction 

Recently, causal viscous hydrodynamics for relativistic heavy ion collisions has been expe- 
riencing rapid development. Several groups independently developed 2+1-dimensional viscous 
hydrodynamic codes that implement longitudinal boost invariance but allow for arbitrary dy- 
namics in the two dimensional plane transverse to the beam So far, most of 
these works have focussed on studying the effects caused by shear viscosity and on developing 
strategies for constraining the QGP shear viscosity to entropy density ratio 77/s from experimen- 
tal data maaaaSaSlil]. ^ was f° unc l tnat elliptic flow vi is very sensitive to tj/s and 
that, due to the rapid expansion of the fireballs created in heavy-ion collisions, even the mini- 
mal KSS bound n/s = l/4n d leads to a strong suppression of compared to the ideal fluid 
case llil 12, Ij, Ifll bl |6|] . A first attempt by Luzum and Romatschke [5] to extract the QGP viscos- 
ity from experimental elliptic flow data, using viscous hydrodynamics simulations, indicate that 
~^\qgp < 5 x j^. To obtain a more precise value requires additional theoretical effort in at least the 
following four directions (see [11] for references): (1) connecting viscous hydrodynamics to a 
hadron cascade to properly account for effects from the highly viscous hadronic stage; (2) includ- 
ing the effects from bulk viscosity; (3) employing a more realistic equation of state (EOS) that 
uses the latest lattice QCD data above T c matched to a hadron resonance gas in partial chemical 
equilibrium below T c , to properly account for chemical freeze-out at r c h e m — 165 - 170 Me V; 
and (4) a better treatment of the initial conditions that not only aims to eliminate presently large 
uncertainties in the initial fireball eccentricity but also properly accounts for pre-equilibrium 
transverse flow and fluctuations in the initial fireball deformation and orientation. 

In this contribution we focus on point (2) and study effects from bulk viscosity. A longer 
account of this work can be found in Ref. B12I1 . In 111 ill , where we reported first results and 
to which we refer the reader for additional details, we constructed and used a function for the 
specific bulk viscosity £/s(T) that interpolates between a "minimal" value well above T c , based 
on strong-coupling calculations using the AdS/CFT correspondence [ 13], to a zero value in the 
hadron resonance gas well below T c , using a Gaussian function between these limits that peaks at 
T c . (For a discussion of uncertainties in (/s below T c see B1411 .) We call this function "minimal 
bulk viscosity". To study larger bulk viscosities, we made comparison runs where this entire 
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function was multiplied with a coefficient C > 1. For twice the minimal bulk viscosity we 
found that the viscous x>2 suppression increases from 20% (for a fluid that has only minimal 
shear viscosity - = ^-) to 30% (for a fluid that additionally features bulk viscosity at twice the 
"minimal" level, C = 2). For a given measured value of V2, accounting for such a 50% relative 
increase in the viscous V2 suppression translates into a reduction of the extracted shear viscosity 
rjl s by roughly a factor |. Consequently, bulk viscous effects cannot be ignored when extracting 
j]/s from experimental data. 

The analysis ifTHl suggests that the main uncertainty stems from insufficient knowledge of 
the peak value of £/s near the phase transition^ However, this is only part of the story. Addi- 
tional complications arise from the fact that the expected peak in (/s near T c is due to rapidly 
growing correlation lengths, associated with a "critical slowing down" of microscopic relaxation 
processes near the phase transition. This probably leads near T c to a much larger relaxation time 
Tn for the bulk viscous pressure than commonly used for the shear viscous pressure (which near 
T c is of the order of 0.5 fm/c, as obtained from kinetic theory JjJ, [H H3] or AdS/CFT JU]). 
Large relaxation times can lead to strong memory effects, i.e. strong sensitivity to the initial 



conditions for the bulk viscous pressure. This is what we discuss here 11211 . 



2. Setup 

We use the code VISH2+1 J2,[3tl to solve (2+ l)-dimensional transport equations for the energy 
momentum tensor T mn and the bulk viscous pressure II: 

d m T" m = , T mn = eu m u n - (p + H)A mn , ( 1 ) 

1 ^ . (rn k 



Dn = (j\ + 0-u)--l\^—d k \-^-u K \. (2) 

Tn 2 r n \(T ) 

Here^m, n denote components in (r, x, y, j] s ) coordinates, with covariant derivative d m (for details 
see iQ), D = u m d m and V ffl = A ml di (where A" 1 " = g" m -u m u n is the projector transverse to the 



flow vector u'") are the time derivative and spatial gradient in the local comoving frame, £ is bulk 
viscosity, and Tn is the corresponding relaxation time. For -(T) we use the phenomenological 
construction described in [11]. For rn we consider three choices: the constant values rn = 0.5 
and 5 fm/c, and the temperature dependent function Tn(T) = max[f • ^(T),0.1 fm/c] with f — 
120 fm/c. The last choice implements phenomenologically the concept of critical slowing down; 
it yields r n * 0.6 fm/c at T = 350 Me V and r n * 5 fm/c at T c . 

To study memory effects, we explore two different initializations for the bulk viscous pres- 
sure: (a) Navier-Stokes (N-S) initialization, II(to) = -(d-u, and (b) zero initialization, II(to) = 0. 
We use To = 0.6 fm/c. For all other inputs we make standard choices as discussed in Refs. llHO] 
and listed in the figure below. 

3. Bulk viscosity effects: uncertainties from relaxation time and bulk pressure initialization 

The left panel of Fig. Q] shows the differential elliptic flow V2(pr) of directly emitted pions 
(without resonance decays) for non-central Au+Au collisions at b — 7 fm, calculated from ideal 



'Currently, theoretical uncertainties for the peak value of | near T c are very large. Extraction from lattice QCD 
simulations gives a peak value around 0.7 flUl . This is more than 10 time larger than the string theory prediction based 
on holographic models 11611 . A critical discussion of the lattice QCD based extraction can be found in fl7ll . 
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Figure 1: (color online) Left: Differential pion elliptic flow v'2(f>r) from ideal and viscous hydrodynamics, including 
only bulk viscosity. Right: Time evolution of the bulk pressure (fl) averaged over the transverse plane (weighted by 
the energy density) from viscous hydrodynamics. Different curves correspond to different initializations and relaxation 
times, as indicated (see text for discussion). 

hydrodynamics and minimally bulk viscous hydrodynamics with identical initial and final con- 
ditions. The different lines from viscous hydrodynamics correspond to different relaxation times 
rn and different initializations II(to). One sees that these different inputs can lead to large un- 
certainties for the bulk viscous V2 suppression. For minimal bulk viscosity, the V2 suppression at 
Pt = 0.5 GeV ranges from ~ 2% to ~ 10% compared to ideal hydrodynamics (blue dashed line 
in the left panel). 

For the shorter relaxation time, rn = 0.5 fmc/c, the bulk viscous v-i suppression is insensitive 
to the initialization of n, and both N-S and zero initializations show ~ 8% \'i suppression relative 
to ideal fluids. The reason behind this becomes apparent in the right panel showing the time 
evolution of the average bulk pressure (II). For short relaxation times, (II) quickly loses all 
memory of its initial value, relaxing in both cases to the same trajectory after about 1-2 fm/c 
(i.e. after a few times Tn). This is similar to what we found for shear viscosity where the 
microscopic relaxation times are better known and short (t„(T c ) ^ 0.2 - 0.5 fm/c) and where the 
shear pressure tensor jt'"" therefore also loses memory of its initialization after about 1 fm/c @1. 

This changes if one accounts for the critical slowing down of the evolution of II near T c . If 
one simply multiplies the constant relaxation time by a factor 10, setting Tn = 5 fm/c, one obtains 
the dotted and solid magenta lines in Fig.Q] Now the bulk viscous v-i suppression relative to the 
ideal fluid becomes very sensitive to the initialization of the bulk viscous pressure: For zero 
initialization II(to) = 0, the viscous vi suppression is very small (only ~ 2% at pj = 0.5 GeV/c). 
The right panel shows that in this case the magnitude of the (transversally averaged) bulk pressure 
evolves very slowly and always stays small, leading to almost ideal fluid evolution. On the other 
hand, if II is initialized with its Navier-Stokes value, which initially is large due to the strong 
longitudinal expansion, it decays initially more slowly than for the shorter relaxation time. Its 
braking effect on the flow evolution is therefore bigger, resulting in much stronger suppression 
of V2 than for zero initialization, at ~ 10% slightly exceeding even the viscous V2 suppression 
seen for the tenfold shorter relaxation time. 

The "critical slowing down" scenario with temperature-dependent Tn(T) (black lines) inter- 
polates between the short and long relaxation times. As for the fixed larger value Tn = 5 fm/c, V2 
depends sensitively on the initialization of IT, but for N-S initialization the viscous v>2 suppression 
is somewhat smaller than for both short and long fixed relaxation times. The reasons for this are 
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subtle since now, at early times, the bulk viscous pressure II evolves on very different time scales 
in the dense core and dilute edge regions of the fireball. As a result, for N-S initialization the 
average value (II) is smaller in magnitude than for both short and long fixed Tn, throughout the 
fireball evolution (right panel, black lines). 

4. Conclusions 

Relaxation times and initial values for the dissipative flows (bulk and shear pressure) are 
required inputs in viscous hydrodynamic calculations, in addition to the transport coefficients 
and the EOS. Near T c , the bulk viscosity f can exceed the shear viscosity 77 of the strongly 
interacting matter. If the relaxation time Tn for the bulk viscous pressure fl is short, it quickly 
loses memory of its initial valaue, but the relatively large peak value of £/s near T c can lead to 
a significant viscous suppression of the elliptic flow V2, competing with shear viscous effects. 
If Tn grows rapidly near T c , due to critical slowing down, the bulk viscous suppression effects 
on V2 depend crucially on the initial value of II: If II is zero initially, bulk viscous effects on V2 
are almost negligible; if II is initially large, however, as for the case of the N-S initialization, 
it remains relatively large throughout the evolution, suppressing the buildup of elliptic flow at a 
level that again competes with shear viscous effects. Additional research on initial conditions and 
relaxation times for the bulk viscous pressure is therefore necessary for a quantitative extraction 
of T]/s from measured data. 
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